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Abstract 



We have used the CLEO II detector to study the multiplicity of charged 
particles in the decays of B mesons produced at the T(4S) resonance. Using 
a sample of 1.5 x 10^ B meson pairs, we find the mean inclusive charged particle 
multiplicity to be 10.71 ±0.02 ^"[^g for the decay of the pair. This corresponds 
to a mean multiplicity of 5.36±0.0ll|]:J^ for a single B meson. Using the same 
data sample, we have also extracted the mean multiplicities in semileptonic 
and nonleptonic decays. We measure a mean of 7.82 it O.OSI'Iq'^q charged 
particles per BB decay when both mesons decay semileptonically. When 
neither B meson decays semileptonically, we measure a mean charged particle 
multiplicity of 11.62 ± OM^Ht Per BB pair. 
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I. INTRODUCTION 



Measurements of the charged track muthphcity distribution in B meson decay are used 
to constrain unmeasured or poorly measured branching fractions in Monte Carlo simulations 
so that generated event samples more closely represent actual data. The CLEO Monte Carlo 
parameterization of B meson decays has been tuned to agree with our measurements and our 
model is used by other experimental groups [||]. Charged particle multiplicity in heavy meson 
decay has been studied by several groups 0] . In this paper we present a measurement 
of the charged particle multiplicity in inclusive BB decays that is an improvement over 
our previous result [Q. We also present improved measurements of the charged particle 
multiplicities in semileptonic and nonleptonic decays. 

For clarity, we use the term "observed multiplicity" to denote the number of well recon- 
structed charged particle tracks in a given event. We use the term "decay multiplicity" to 
denote the number of e^, fi^, vr^, and that come from the decay of the primary 
B mesons and also from the subsequent decays of any secondary or tertiary particles other 
than neutrons, K^, or 7r°. The decay multiplicity excludes any tracks produced through 
interactions with the detector or surrounding material. Not all charged decay products will 
result in reconstructed tracks, and not all observed tracks come from the primary decay, so 
the observed multiplicity may be less than, equal to, or greater than the decay multiplicity 
for a given event. 

II. INCLUSIVE MULTIPLICITY MEASUREMENT 

The CLEO detector is located at the Cornell Electron Storage Ring, a high luminosity 
e'^e~ collider operated at or near the T{4S) resonance. The results presented here are derived 
from a sample of 1.4 fb~^, corresponding to 1.5x10^ B meson pairs, collected with the CLEO 
II detector ||^. Charged particle tracks are measured by cylindrical wire drift chambers inside 
a 1.5 T superconducting solenoid. A Csl crystal calorimeter is also inside the magnet, and 
energy deposition information from both the calorimeter and the drift chamber is used for 
particle identification. Muon counters are layered in the steel yoke surrounding the coil. 

To obtain a clean sample of candidate BB events, we select hadronic events by requiring 
that an event have three or more reconstructed tracks, energy deposition in the calorimeter 
greater than 15% of the center of mass energy and an event vertex consistent with the 
interaction region. For additional background suppression, the total reconstructed event 
energy, including charged and neutral particles, is required to be between 4 GeV and 12 GeV, 
and the total reconstructed vector momentum of the event is required to have a magnitude 
less than 3 GeV/c. This hadronic event sample contains events from both BB and continuum 
processes such as qq and r+r^ production. We remove the continuum contribution by 
rescaling and subtracting the observed multiplicity distribution of a separate 0.7 fb~^ data 
sample collected 65 MeV/c^ below the T(4S') resonance. 

To be counted in our observed multiplicity, drift chamber tracks are required to be well 
reconstructed and consistent with having originated from the event vertex. Tracks must 
not be within 25.8 degrees of the e~^e~ beam axis. Once the event selection, continuum 
subtraction, and track selection are completed, we count the selected tracks in each event to 
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obtain the observed charged track muhiphcity distribution in Fig. |T|. There are events with 
fewer than three selected tracks because not all reconstructed tracks pass the track selection 
criteria. 
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FIG. 1. Observed charged track multiplicity for BB events. 

To obtain the true decay multiplicity distribution from the observed multiplicity distribu- 
tion, we must account for detector effects. If the number of events with observed multiplicity 
j is Oj and the number of events with decay multiplicity i is Di, these quantities are related 
by Eqn. |1]. 



E 



i=2A,. 



where eij is the probability that an event with decay multiplicity i will be reconstructed with 
observed multiplicity j. We have assumed that charge is conserved, so the index i can only 
take even values. In principle, there can be events with decay multiplicity of zero, where two 
neutral B mesons decay to all neutral final states. However, we do not include zero decay 
multiplicity events in our analysis both because of the very low branching ratio for such 
events and also because our event selection criteria make detection of such events extremely 
unlikely. The upper bound in Eqn. |1| is, in principle, the maximum decay multiplicity in a 
BB event, which is not known. We vary the maximum decay multiplicity in our analysis as 
described below. The fact that the sum of and the sum of Oj are both equal to the total 
number of events is used to constrain the values of the Di, as expressed in Eqn. ^. 



E A = Eo, 



(2) 



=2,4,. 



i=0 
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where the upper bound m is the maximum value of our observed multipicity which is 20. 

The coefficients e^j in Eqn. |I] are obtained from Monte Carlo simulation and depend 
primarily on the detector's track finding efficiency and also on the probability of producing 
extra charged particles that pass the track selection cuts and that are counted. While these 
coefficients depend on accurate simulation of detector response and processes such as photon 
conversion and decays in flight, they do not depend significantly on the exact tuning of the 
branching fractions or the decay multiplicity distribution in the simulation. 

The parameters Di in Eqn. |l] are determined by a fit- This fit unfolds the detec- 
tor effects to give the decay multiplicity distribution of events that pass the event selection 
criteria. These selection criteria are biased against very low decay multiplicity events, par- 
ticularly because of the requirement of three or more reconstructed tracks. We remove the 
event selection bias using Monte Carlo simulation to determine the probability for events 
of a given decay multiplicity to pass the event selection cuts. After unfolding detector and 
reconstruction effects, we obtain the decay multiplicity distribution in Fig. ^ The error bars 
represent the statistical uncertainty in both the fit and the event bias correction, but 
do not include systematic errors. The dashed lines represent the high multiplicity and low 
multiplicity statistical fluctuations in the fit. The large error bar on the i = 2 point is due 
to the event bias correction. 
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FIG. 2. Unfolded charged track decay multiplicity for BB events. The dashed lines represent 
the high multiplicity and low multiplicity fluctuations in the fit. 



From the distribution in Fig. 0, we obtain a mean of 10.71 ± 0.02 charged particles for 
inclusive BB decay, where the error is statistical only. The points with i > 20 show no 
evidence for events with such high multiplicities. We have tested our fitting procedure by 
varying the maximum decay multiplicity included in the fit between 20 and 28. The fit is 
stable and the unfolded mean decay multiplicity does not change significantly when decay 
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multiplicities of 22 and higher are included or excluded from the fit. 

The most important systematic effect in this analysis is the accuracy of modeling the 
detector's track finding efficiency. Our studies indicate that the overall efficiency for finding 
a single track is known to within ±1% for tracks with momentum greater than 250 MeV/c 
and with decreasing accuracy as track momentum decreases. The uncertainty in single track 
finding efficiency gives a 11qI% uncertainty in the measured mean decay multiplicity. Removal 
of the event selection bias shifts the measured mean decay multiplicity by +0.7%, which we 
also take as the systematic uncertainty for event selection bias. Because the track finding 
efficiency depends on the momentum of the track, we account for the uncertainty due to this 
dependence. This analysis uses all tracks without regard to momentum. When we add a 
track selection cut requiring a reconstructed momentum of at least 150 MeV/c, we observe a 
shift in the decay multiplicity of —0.6%. We assign an additional ±0.6% uncertainty due to 
this momentum dependence. Charged pions produced by the decays of Ks will have lower 
track finding efficiency than average charged pions, so our result depends on the rate of Ks 
production in our simulation. Based on our studies of inclusive Ks production in B meson 
decay, we assign a systematic uncertainty of Iq 3% from this source. Interactions of neutral 
and charged decay products with the detector material produce additional charged tracks, 
some of which satisfy the track selection criteria. We study the effect of these extra tracks 
by varying the rates of photon conversion and hadronic interactions in our Monte Carlo 
sample. Misidentification of these extra tracks contributes an additional ±0.3% systematic 
uncertainty. Extra tracks also come from the decay of particles in the detector volume, 
adding another ±0.1% uncertainty. Additional uncertainty comes from contamination by 
non-BB events, most notably from beam-gas interactions. We study the effect of non-BB 
events by varying the size of the continuum subtraction by 1%, which is the uncertainty in 
our measured luminosity, and we find a ±0.2% uncertainty in our mean decay multiplicity. 

We have looked for other potential systematic effects by varying our track selection cuts 
but do not observe any significant change in mean decay multiplicity. When all systematic 
errors are added in quadrature, we obtain a total systematic uncertainty of 1l4%, which gives 
a final result of 10.71 ± 0.02lgj5 for the mean inclusive charged particle decay multiplicity 
in the decay of a BB pair. Systematic errors are summarized in Table |. 

The results presented here are a significant improvement over the previous CLEO re- 
sult [0. The earlier analysis found a mean multiplicity of 11.5 ± 0.2 ± 0.4. In addition to 
using a much larger sample of BB events, the present analysis uses a substantially different 
method. The former analysis treated the coefficients etj in Eqn. as a matrix and used 
matrix inversion to solve for the Di. Such matrix inversion is very sensitive to singularities 
and can produce unstable results. In the current analysis, we do not attempt to invert the 
eij matrix, but instead we perform a more stable fif fo^' fhe parameters Di. 



III. SEMILEPTONIC AND NONLEPTONIC MULTIPLICITY MEASUREMENT 

We have further analyzed the same data sample to measure multiplicities separately for 
semileptonic and nonleptonic decays. We define semileptonic to include only decays of B 
mesons into an electron or muon plus a neutrino and any number of hadrons. All other decays 
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Systematic Error Source 


Effect on Mean Decay Multiplicity 




(Charged Tracks) 


Overall Track Finding Efficiency 


+0.17 
-0.09 


Event Selection Bias 


±0.07 


Low Momentum Tracking Efficiency 


±0.06 


Ks Modeling 


+0.05 
-0.03 


Hadromc interactions and rnoton Conversion 


±0.03 


Contamination by Non BB Events 


±0.02 


Decays in Flight 


±0.01 


All Other Sources 


±0.01 


Total Systematic Error 


+0:11 

-0.15 



TABLE I. Systematic Errors for the Inclusive Decay Multiplicity Measurement for BB Pairs. 



are classified as nonleptonic. Decays involving tau leptons are counted as nonleptonic because 
approximately 65% of taus decay hadronically and such events cannot be distinguished 
easily from purely hadronic B decays. 

We use the same event and track selection criteria as in the inclusive analysis and sort the 
events by the number of leptons identified. We identify electrons by combining information on 
energy deposition in the drift chamber, the shape of the shower observed in the calorimeter, 
and the ratio of the calorimeter energy to the track momentum. Muons are required to 
traverse at least three pion nuclear interaction lengths of material. We also require all lepton 
candidates to have momentum between 1.4 GeV/c and 2.5 GeV/c to suppress false lepton 
identification and secondary leptons. 

Once we have sorted the events by the number of leptons found, we correct for misiden- 
tified and secondary leptons. We move these misidentified events from the sample with one 
identified lepton to the sample of events with no identified lepton by using the observed 
multiplicity distributions of Monte Carlo generated events that are similarly misidentified 
in the reconstruction. The fake and secondary observed multiplicity distributions are scaled 
using our best estimates of the fake and secondary lepton rates and are added to the sample 
with no identified lepton and subtracted from the sample with one identified lepton. Un- 
certainties in the rates and observed multiplicity distributions of these events are treated as 
systematic errors. 

We use the same methods as in the inclusive analysis to extract the mean decay mul- 
tiplicity from the corrected observed multiplicity distributions for the samples with zero or 
one identified lepton. For events with no reconstructed leptons, we obtain a mean decay 
multiplicity of 11.02 ± 0.01, and for the sample with one reconstructed lepton, we find a 
mean of 9.28 ± 0.02 charged tracks, where the errors are statistical only. 

Not all semileptonic decays will produce a detected lepton because not all electrons and 
muons will enter the fiducial tracking volume of the detector and pass the track selection 
and lepton identification criteria. To obtain the true decay multiplicities for semileptonic 
and nonleptonic BB decays, we use simulation to unfold the migration between the number 
of leptons generated by the decay of the B mesons and the number of leptons actually 
identified. If the migration probability for a given event were independent of the event's 
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decay multiplicity, we could write: 

= a'Mn-n + b'M,_,, (3) 

where the M* are the mean decay multiplicities of the sample with i identified leptons, 
Mn-n and Ms_s are the mean decay multiplicities of events where both B mesons decay 
nonleptonically and semileptonically, respectively, and a* and 6* are migration coefficients 
determined from Monte Carlo simulation. These multiplicities include the lepton tracks. 

There are many factors that affect the number of leptons that will be reconstructed for a 
given event. Lepton identification depends strongly on particle momentum not only because 
of the explicit momentum requirements but also because a particle's momentum determines 
the probability that it will penetrate the muon chambers sufficiently to be classified as a muon 
candidate. Because of phase space limitations, events with one or more high momentum 
(above 1.4 GeV/c) particles will tend to have fewer charged tracks than events without a 
high momentum particle, so events with an identified lepton tend to have lower multiplicities 
than those where the lepton was not identified. 

In addition to assuming that migration probability is independent of decay multiplicity, 
Eqn. ^ also assumes that the mean decay multiplicity of events where only one B meson 
decays semileptonically is simply the average of Mn-n and Mg^g. This assumption would be 
valid if the mean true decay multiplicities oi B~ and B^ mesons were equal, but this may 
not be the case. To account for the possibility of unequal B~ and B^ multiplicities and the 
effect of multiplicity dependent migration, we can introduce another term into Eqn. |^: 

M' + A' = dMn-n + h'Ms_s. (4) 

The correction terms A* and migration coefficients are determined from Monte Carlo 
simulation, and we can solve Eqn. ^ for M„_„ and M^.^. We obtain mean true decay 
multiplicities of 11.62 ± 0.04 charged particles for events where neither B meson decays 
semileptonically and 7.82 ± 0.05 charged particles for events where both B mesons decay 
semileptonically, where the errors are statistical and include analytic propagation of the 
uncertainty in each of the parameters in Eqn. ^. 

All of the systematic uncertainties that affect the inclusive analysis are also present in 
these semileptonic and nonleptonic results. This li;4% uncertainty in unfolding the decay 
multiplicity from the observed multiplicity is the largest systematic effect for these results. 
The next largest uncertainty comes from lepton identification, including inefficiencies, false 
particle identification, and acceptance of secondary leptons not directly produced in the 
B meson decay. We have studied these effects by varying all of our lepton identification 
criteria. We have also varied the rates of fake and secondary lepton identification directly 
in simulation. All of the lepton identification systematics combine to give a ±0.06 charged 
track uncertainty in the nonleptonic result and a corresponding uncertainty of ±0.11 charged 
tracks in the mean decay multiplicity of semileptonic decays. Our result is also sensitive to the 
modeling of semileptonic decays in our simulation, including the inclusive rate of semileptonic 
decays and the exclusive branching fractions of semileptonic decays that produce the various 
resonances of the D meson. We have varied both the inclusive and exclusive branching 
fractions within their uncertainties and find ±0.05 uncertainty in the nonleptonic mean 
decay multiplicity and ±0.08 uncertainty for the semileptonic result. Uncertainty in the 
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Systematic Error Source 


Nonleptonic Events 


Semileptonic Events 




(Charged Tracks) 


(Charged Tracks) 


Inclusive Unfolding 


-0.16 


+0.15 
-0.11 


Lepton Identification 


±0.06 


±0.11 


Branching Fractions 


±0.05 


±0.08 


Fake and Secondary Leptons 


±0.01 


±0.06 


All Other Sources 


< ±0.01 


±0.01 


Total Uncertainty 


+0.24 
-0.18 


+0.21 
-0.19 



TABLE II. Systematic Error for the Charged Track Mutliplicity in Nonleptonic and Semilep- 
tonic Events. 



observed multiplicity of events with fake and secondary leptons contributes ±0.06 charged 
tracks uncertainty to the semileptonic result. The systematic uncertainties are summarized 
in Table 0. 

Combining all errors, we find mean multiplicities of 11.62 ±0.04^gj8 charged tracks when 
both B mesons decay nonleptonically and 7.82±0.05loj9 charged tracks when both B mesons 
decay semileptonically. As a further check on our systematic errors, we have repeated this 
analysis by using only electron identification or only muon identification. In both cases, the 
results are statistically consistent with the means stated above and with each other. Our 
present results are a significant improvement over the previous CLEO results 0, in which 
mean multiplicities of 12.6 ± 0.4 ± 0.4 and 8.2 ± 0.7 ± 0.4 were obtained for events where 
both B mesons decay nonleptonically and semileptonically, respectively. 
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